Alignment and scaling of large-scale fluctuations in the solar wind 
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We investigate the dependence of solar wind fluctuations measured by the Wind spacecraft on 
scale and on the degree of alignment between oppositely directed Elsasser fields. This alignment 
controls the strength of the non-linear interactions and, therefore, the turbulence. We find that at 
scales larger than the outer scale of the turbulence the Elsasser fluctuations become on average more 
anti-aligned as the outer scale is approached from above. Conditioning structure functions using the 
alignment angle reveals turbulent scaling of unaligned fluctuations at scales previously believed to 
lie outside the turbulent cascade in the '1// range'. We argue that the 1// range contains a mixture 
of non-interacting anti-aligned population of Alfven waves and magnetic force-free structures plus 
a subdominant population of unaligned cascading turbulent fluctuations. 



Introduction. The solar wind is a hot, tenuous plasma 
that flows away from the Sun at supersonic speeds. Tur- 
bulence transports energy from the driving 'outer' scale 
to smaller scales via non-linear magnctohydrodynamic 
(MHD) interactions of magnetic B and velocity V fields, 
until kinetic effects and dissipation become important 
close to the ion gyroscale. In fast solar wind ( [ | > 600 
km/s), a '1//' scaling of magnetic-field power spectra is 
observed at low spacecraft frequencies, / HH3. Slowly 
evolving structures are advected at supersonic speeds 
past spacecraft, so the observed spacecraft frequency of 
a fluctuation is proportional to its characteristic wave 
number (scale) k The energy spectrum in the 1/f 
range is, therefore, expected to scale as E{k) oc fc _1 . 
A steeper spectrum close to fc -5 / 3 associated with tur- 
bulence is observed at higher spacecraft frequencies in 
the 'inertial range' and there is a spectral break between 
the two regimes 0, H, H ! at 1 AU this typically oc- 
curs at / ~ 10 -3 Hz. Studies have shown 0, y|, @- 
IIo| that the power spectral density of fluctuations in 
the low-frequency band decreases with distance from the 
Sun as i? -3 , consistent with these scales containing non- 
interacting Alfven waves ll| , and thus these large-scale 
fluctuations are thought to have originated at the Sun 
and travelled outwards with relatively little in- situ modi- 
fication. In this Letter, we argue that this interpretation 
is incomplete as nonlinear interactions occur at larger 
scales than previously thought (see also (l2|). 

Recently the concept of scale-dependent alignment has 
become prominent in theoretical and numerical studies of 
MHD turbulence [IM!]. Scale-dependent alignment is 



the tendency for the angle between fluctuations of B and 
V in the plane perpendicular to the mean magnetic field 
Bo to decrease with increasing k. Attempts to measure 
alignment in the inertial range of the solar wind produce 
no evidence of scaling [lj| [2(J, although the ability to 



measure the scale dependence is limited by instrument 
noise characteristics. These studies do, however, find a 
scaling of the alignment in the 1/f range, which is unex- 
pected given the previous interpretations of these large- 
scale fluctuations as non-turbulent. 

Here we study the alignment of Alfvenic fl uctuati ons 
in the 1/f range. 
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We use Elsasser variables 
characterize the Alfven waves that travel sunward and 
anti-sunward in the plasma frame. We define the angle 
cf> between fluctuations in the Elsasser fields in the plane 
perpendicular to Bo as the alignment angle. This angle 
is geometrically related to the alignment angle between 
fluctuations <5B and 5V but is not completely determined 
by it. The alignment angle (f> is important because it con- 
trols the strength of the non-linear interaction [24| and, 
as we are about to see, allows one to sort the large-scale 
fluctuations into steep-scaling 'turbulent' and shallow- 
scaling 'non-turbulent' populations. 

Data. We use Wind spacecraft observations of solar 
wind magnetic field B, velocity V, and proton number 
density n p at cadence 5t = 3 s made by the MFI and 3DP 
instruments during a 6-day-long fast stream interval ob- 
served between days 14 and 20 of 2008. The average so- 
lar wind conditions were: \V\ = 660 km/s, \B\ = 4.4 nT, 
n p = 2.4 cm~ 3 , Alfven speed Va = 62 km/s, and the ra- 
tio of thermal to magnetic pressure for protons (3 P = 1.2. 
Similar fast streams recurred five times in succession due 
to a long-lived, low-latitude coronal hole. A further three 
of these fast streams were also analyzed and provide 
quantitatively similar results to those shown here; one 
stream was excluded because it coincided with a large 
data gap. The same analysis performed on Ulysses space- 
craft data when in fast wind over the poles of the Sun also 
shows qualitatively the same results as described below. 

We use the Alfven-normalized magnetic field B = 
B/ yJAirm p n p . Elsasser variables, = V ± B, are re- 
defined so that Z + are anti-sunward and Z~ sunward 
propagating fluctuations in the plasma frame. We are 
interested in alignment, so we use only the projection of 
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FIG. 1. Structure functions of Wind data during a 7-day- 
long fast stream. Tg is the estimated travel time of the solar 
wind from the Sun to the spacecraft, Ta is the time it takes 
for the largest scale over which Alfven waves can interact to 
be advected past the spacecraft, and To is the time scale 
corresponding to the outer scale of the turbulence, defined as 
the spectral break in the B structure functions. The bottom 
panel shows normalized cross helicity a c and residual energy 
oy calculated from this data. 



the fluctuating fields on to the plane perpendicular to the 
local mean magnetic field Bq at a time scale r: 



6t 



t'=t + T 



B (t,r)=- ]T B(f), 



(1) 



5x(t,r)=X(t)-X(t + r), (2) 
5x±(t,T) =Sx(t,r) - (Sx(t,T)-B Q (t,T)) B (t,r), (3) 



where * denotes unit vectors and X can be B, V, Z + 
or Z . In the plots presented below, the time scale r 
is converted into a frequency in the spacecraft frame to 
facilitate comparison with Fourier spectra: / = 1/r. A 
logarithmically spaced range of time scales 6 s < r < 
2 x 10 5 s is used to investigate the inertial and 1// ranges 
of the fast solar wind. 

Structure functions. In Fig. [TJ we show the second- 
order structure functions of all four vector fields perpen- 
dicular to the magnetic field: 



S 2 (<5x,r) = 1 £ |<5x ± (t,r)| 2 = (|<Jx ± (t, r)| 2 
i=ti 



where N is the number of samples in the time period t\ < 
t < ti. The scaling exponent of the structure functions 
a, where S^Jx, t) oc r~ a oc is related to the Fourier 
spectral index 7 by 7 = a — 1 [SB], |2(| . 

The vertical lines in Fig. [1] show important time scales 
for this period of solar wind. T s = 1 AU/|F| = 2.3xl0 5 s 



FIG. 2. The mean (black circles) and the probability distribu- 
tion (color scale) of <f> at different scales. The angle 9 between 
(5vj_ and Sh± is shown as white squares for comparison. 



is the approximate time the solar wind has taken to travel 
from the Sun to the Wind spacecraft. To — 2 x 10 3 s 
is the approximate time scale associated with the outer 
scale, defined as the scale at which S% (5h) rolls over from 
flat in Fig. [T] (a = corresponding to the spectral in- 
dex 7 = — 1, the 1/f range) to an inertial range scaling 
(a w -2/3). T A = L/\V\ = 5.4 x 10 4 s is the approximate 
time scale associated with the advection past the space- 
craft of the largest separation L two counter-propagating 
Alfven waves can have and still meet one another in the 
time the solar wind has taken to propagate from the Sun 
to the spacecraft. To calculate Ta we estimate the depen- 
dence on heliocentric distance R as follows: \B\ cx i? -15 
and pi cx R~ 2 and thus Va cx R~ - 5 and solve for the dis- 
tance from the Sun L that the slower of the two Alfven 
waves (\V\ + Va and \V\ — Va) must start so that the 
faster wave just meets it at 1 AU. Thus the spacecraft 
frequencies between / ~ I/T5 and / ~ 1/Ta represent 
spatial structure between different source regions in the 
corona, since they cannot have interacted during tran- 
sit from the Sun. The range of frequencies between 1/Ta 
and 1 /To contains fluctuations that may have interacted; 
on these scales, all structure functions are relatively flat, 
with S2(Sbj_, t) and S2{5z\, r) oc f°. Frequencies higher 
than 1 /To show all variables with scaling typical of tur- 
bulence in the fast solar wind: 5 2 , r) o c Sn (Szj_ , r) oc 



(4) S 2 (<5zI,t)oc/- 2 / 3 ,S 2 (5v,t)cx/- 



Alignment angle. In order to investigate correlations 
between Elsasser fluctuations, we calculate the local 
scale-dependent <f>: 



<t>(t,T) 



arccos 



6z+(t,T) ■ SzZ(t,T) 
\Sz+(t,r)\ \Sz2(t,r)\ 



(5) 



In Fig. [21 we show the mean (in black) and the probability 
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FIG. 3. Conditioned structure functions of Elsasser variables 
from observations of fast solar wind. Anti-aligned (circles) 
and un-aligned (squares) structure functions of 8z~\_ (blue) 
and Sz^_ (green) are plotted against frequency. Un-aligned 
fluctuations begin scaling at / ~ 1/Ta and have similar gra- 
dients to those found in the inertial range, indicating that 
they may be interacting non- linearly in the 1// range. 

distribution of </> at each scale. The distribution of 4> has 
been discretized using 10° wide bins. At all scales, the 
distribution of (j> covers the full range of possible values 
and is peaked at 180°. The mean is not the most probable 
value at any scale and is strongly dependent on the tail 
of the distribution that extends towards 0°. The mean 
values of 8, the angle between 5v± and <5bj_, calculated 
in a similar manner to </>, are also shown (in white). 

The frequencies / ~ 1/Ta and / ~ l/To both coin- 
cide with marked changes in behavior of the distribution 
and mean of <f> in Fig. [5J The mean alignment angle, 
(</>), increases between / ~ 1/Ta and / ~ l/To, but 
then rolls over and decreases at higher frequencies, sim- 
ilar to previous observations of 9 fljil . l2oT ] . The peak in 
the distribution at 180° grows as frequency increases in 
the range 1 /Ta < / < 1 /To but then flattens and begins 
to decrease where instrument noise becomes important, 
as discussed below. 

Relating alignment and spectral scaling. To investigate 
whether </> has any effect on the turbulence and what may 
be causing the change in <j) with scale, we use structure 
functions conditioned on <f>. Eq. §4§ is modified to average 
over instances when <fi is within a fixed range: 

0(t,r)=0 2 

S 2 (Sz ± ,r\4>i<<t><h) = ^ E l 5z t(*. T )f. ( 6 ) 

0(t,r)=0i 

where N is the number of points with <pi < <p < <f> 2 . 

The structure functions calculated according to Eq. (|6]) 
are shown in Fig. [3] for the anti-aligned and perpendicu- 



l/T A < f < l/To l/To < / < 5 x 1(T 3 

80<<ft< 90 170 <(j>< 180 80 <0< 90 170 <(/>< 180 
8z\ -0.32 ±0.02 -0.07 ±0.02 -0.49 ± 0.06 -0.36 ± 0.02 
8z~ -0.65 ± 0.06 -0.22 ±0.03 -0.60 ± 0.05 -0.43 ± 0.02 



TABLE I. Scaling of the structure functions S2(Sz ± ,r\(f>) oc 
f a of the Elsasser fluctuations in two frequency ranges cov- 
ering the 1/f range and the inertial range. The values are 
calculated from a linear least squares fit of a straight line to 
the structure functions on a log- log plot. 

larly aligned fluctuations: S 2 (Sz ± , t|170° < <f> < 180°) 
(filled symbols) and S 2 (Sz ± ,t\80° < <j) < 90°) (open 
symbols), respectively This allows us to investigate the 
effect that the alignment angle has on the turbulence by 
separating differently aligned fluctuations and observing 
their scaling in the 1/f and inertial ranges. 

Between / ~ 1/Ta and / ~ l/To, the anti-aligned 
anti-sunward <5zJ structure functions scale with a = 
— 0.07± 0.02, giving a spectral index close to —1. These 
fluctuations are the most common (Fig. [2| and contain 
the most power and hence dominate the bulk average 
structure functions in Fig. Q] - and presumably all pre- 
viously reported, unconditioned structure functions and 
spectra in this frequency range of fast solar wind. Per- 
pendicularly aligned 5z\_, however, have a steeper scaling 
a = —0.32 ± 0.02, corresponding to a spectral index only 
slightly shallower than that at higher frequencies in the 
inertial range, a = —0.49 ± 0.06. 

The structure functions of perpendicularly aligned sun- 
ward fluctuations 6z^_ are steep from / ~ 1/Ta until the 
instrument noise floor (the solid red line; see discussion 
below) is reached, with a = —0.65 ± 0.06 in the 1/f 
range and a = —0.60 ± 0.05 in the inertial range, giving 
a spectral index close to —5/3 in both frequency ranges. 
The anti-aligned <5zJ structure functions are flatter in 
the range 1/Ta < f < l/To than in the inertial range, 
with a = -0.22 ±0.03 and a = -0.43 ±0.02 respectively 
These scalings arc summarized in Table |U 

Accuracy of measurements. Measurement noise is a 
potential concern in this analysis. The 3DP instrument 
is known to have noise in the high-cadence moments [l9T | 
with observations often appearing discretized. By differ- 
encing the raw velocity data and finding the most com- 
mon value we estimate the noise amplitude during the 
periods we analyze to be approximately equivalent to a 
2 km/s uncertainty in each component of V. A standard 
error analysis on Equation (|4]) leads to the frequency de- 
pendent noise represented by the red line in Fig. [3] This 
noise affects the structure functions of perpendicularly 
aligned sunward fluctuations the most, with the signal 
to noise ratio of un-aligned <5zJ structure functions be- 
coming sig nificant (~ 2) at / ~ 2 x 10~ 3 Hz; this is our 
estimate of the frequency at which noise begins to ren- 
der our results unreliable. We have, therefore, restricted 
the fitting of the structure functions scaling in the in- 
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ertial range (Table [XJ) to the lowest frequency decade, 
1/T < f < 5 x 10~ 3 Hz. This minimizes the effect of 
the noise on the scaling and we do not draw conclusions 
about the inertial range. 

The noise will also affect the measurement of cj) since 
it uses the values of 6z~J_ and 6z~^ , which in turn contain 
V observations. The roll-over of (cj)) at / ~ 1/Tp in 
Fig. [5] occurs at a frequency a factor of 4 lower than the 
frequency /~2xlO _3 Hz, at which the signal to noise 
ratio of the weakest structure functions reaches a value of 
2, and so cannot solely be attributed to noise. The strong 
decrease in (cj)} and the flattening of the distribution at 
spacecraft frequencies above 2 x 1CP 3 Hz is, however, 
likely to be caused by the noise. 

Alignment and geometry. Alignment of 5h± and 6v± 
is related to the alignment of 8z\ and <5z]\ Fig. |4] shows 
a typical geometry in the plane perpendicular to Bo in 
the fast solar wind, assuming 5h± > 5v± with only a 
small angle 9 between the vectors, as suggested by Fig. [I] 
and [2j This results in 6z~^ > <JzJ with a large angle 
4> between them (anti-alignment). This simple geom- 
etry can be expressed in terms of the scale-dependent 
dimensionless parameters (29| normalized cross helicity 
cj c = (\6z\\ 2 - |<5z-| 2 )/(|(5z+| 2 + |5zX| 2 ) and normalized 



residual energy oy = ( 1 6v± | 2 — | <5b^ 2 ) / ( | 6v j 



|<5bj 



cos(</>) = 
cos(0) = 



oy 



6z^_ ■ 5zj_ 

\S*1\\6*1\ 
<5vj_ ■ 6b± oy 
\6v ± \\6b ± \ = yi - ^ 



(7) 
(8) 



The geometry in Fig. [4] is fixed by setting any two of 
9, (j>, a c and oy. Therefore, statements about alignment 
are also statements about normalized cross helicity and 
residual energy of individual fluctuations, and vice- versa. 
The fluctuations that scale steeply in Fig.[3]are those with 
(j) ~ 90° and therefore oy ~ meaning that 8v\ ~ 5b\. 

While we note that the relations in Eq. ([7J [5]) are only 
strictly true for individual realizations it is interesting to 
look at our results concerning averages from this point of 
view In the range 1/Ta < / < 1/To, o r is negative 
and decreases towards —1 while oy increases towards 1 
(Fig. [T]). So, in accordance with Eq. (UJ, the alignment 
angle cj) tends towards 180° (Fig. [2]). Note that this situ- 
ation is perhaps consistent with the idea that MHD tur- 
bulence would generate negative residual energy In 
the same vein, we conclude from Eq. © that the align- 
ment between i5vj_ and Sb± intensifies in the 1// range, 
as indeed seen in the solar wind (the angle 9 is shown in 
white in Fig. [5] and was previously measured in flil I20I]). 

Discussion. We have shown that in the fast solar 
wind, both the distribution and the mean of the an- 
gle between Elsasser fluctuations, <f> are scale depen- 
dent. The probability of Elsasser fluctuations being anti- 
aligned (<f) ~ f 80°) starts increasing at the scale at which 
Alfvcn waves begin to interact, / ~ 1/Ta, and stops at 
the outer scale, / ~ 1/Tq. 




5z + , 



FIG. 4. Fluctuation vectors and alignment angles in the plane 
perpendicular to B. 



In Fig. [3] and the attendant discussion, we showed that 
the I // scaling of spectra (flat scaling of structure func- 
tions) detected in numerous previous observations of fast 
solar wind is dominated by the anti-aligned subset of 
the anti-sunward Elsasser fluctuations (5z~\_). Hidden be- 
neath this energetically dominant sea of non-interacting 
(or weakly interacting) 'non-turbulent' fluctuations are 
the unaligned fluctuations, which exhibit steep spectral 
scalings symptomatic of a nonlinear cascade. We hypoth- 
esize that the steep scaling of unaligned fluctuations is 
caused by the increased in-situ non-linear interaction of 
these fluctuations, since both populations of fluctuations 
have similar cr c and travel time from the Sun. The dif- 
ferent behavior of these populations is reminiscent of the 
difference between slow and fast solar wind streams [35| , 
however the fluctuations that scale steeply (oy = 0) do 
not resemble those characteristic of slow wind (oy < 0). 

It is an interesting question whether the anti-aligned 
'non-turbulent' fluctuations are Alfven waves or magnet- 
ically dominated force-free structures. The pure case of 
the former would require Sz^ 3> <5zJ and so <5bj_ ~ <5v^ 
(oy w 1, |oy| <C I); the pure case of the latter, Sb± ^> 
5v± and so 5z~\_ ~ 5zJ_ (oy rj 1, |oy| <C 1). The mea- 
sured fluctuations appear to be in between these two ex- 
tremes (6b ± > 6v± and <5zJ > 6z ~) a nd can perhaps 
be interpreted as a mixture of them [361 . |37| . Both types 
of fluctuation are slow to decay; this can be thought of 
in terms of conservation of cross helicity (Alfven waves) 
and magnetic helicity (force-free structures, subject to 
the minimum-energy constant-helicity relaxation princi- 
ple 38, [3^). The generation of residual energy at low 
frequencies [34[ could then be interpreted as generation 
(or occurrence and persistence) of force-free structures. 

We conclude that the turbulent cascade in the fast so- 
lar wind starts at larger scales than previously thought, 
although it is restricted to perpendicularly aligned fluctu- 
ations and energetically sub-dominant. Measured scale- 
dependent alignment in the I// range represents the 
change in the fractional populations with scale of tur- 
bulent, non-linearly interacting, perpendicularly aligned 
fluctuations versus non-interacting, anti-aligned fluctu- 
ations. We have identified a new, larger, outer scale 
(/ ~ 1/Ta), which is consistent with an anti-sunward 
Alfven wave requiring only one interaction with an oppo- 
sitely directed wave to launch the perpendicularly aligned 
cascade. It is still uncertain what determines the fre- 
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quency at which the spectral break between the 1// and 
the inertial ranges occurs [H, 0, H3, 12, 
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